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A symmetric dissipative two-state system is asymptoti- 
cally completely delocalized independent of the initial state. 
We show that driving-induced localization at long times can 
take place when both the bias and tunneling coupling energy 
are harmonically modulated. Dynamical symmetry breaking 
on average occurs when the driving frequencies are odd mul- 
tiples of some reference frequency. This effect is universal, 
as it is independent of the dissipative mechanism. Possible 
candidates for an experimental observation are flux tunneling 
in the variable barrier rf SQUID and magnetization tunneling 
in magnetic molecular clusters. 

PACS: 05.30.-d, 05.40+j 



The dissipative two-state system is the simplest model 
allowing the study of dccoherence in a tunneling system 
jl],|]]. The two states may model, e.g., the localized states 
in a tight-binding double well. This model has found 
broad application to tunneling systems in physics and 
chemistry ||||. It has been used, e.g., to study the tun- 
neling of atoms between the tip of an atomic-force micro- 
scope and a surface H,0] , the dynamics of the magnetic 
flux in a superconducting interference device (SQUID) 
P], or the relaxational behavior of tunneling centers in 
submicron Bi-wires |7j . Also the low temperature dynam- 
ics of the magnetization of anisotropic magnetic molecu- 
lar clusters can be described in terms of a two-state sys- 
tem |8|]. Tunneling of magnetization has been observed 
in molecular Mn 12 |J and Fe 8 |1(J nanomagnets. 

Driven dissipative systems exhibit a variety of dynam- 
ics and opened new potentials of application pJj ]. For 
example, the complete suppression of tunneling of a bare 
two-state system (TSS) induced by a suitably chosen 
monochromatic driving field |12| can also persist in the 
presence of dissipation Jl3J — |l 5|] . A variety of different 
noncquilibrium initial preparations, relevant e.g. in elec- 
tron transfer, can equivalently be described in terms of 
a particular time-dependent driving force U&. An ap- 
plied dc-ac field can invert the population of donor and 
acceptor sites 117]]. Time-dependent modulation of the 
bias energy can be achieved, e.g., by an ac electric field 
coupled to the dipole moment of the TSS, or, in a SQUID 
ring device, by a time-dependent external magnetic field 
threading the ring. In contrast, "quadrupole" -like cou- 
plings induce a variation of barrier height and width and 
thus lead to a modulation of the tunneling matrix ele- 
ment (TME). The effects of monochromatic modulation 
of the TME have been studied in Ref. [[18|. Promising 



candidates to study the effects of concerted harmonic 
modulation of the bias and TME are the variable bar- 
rier rf SQUID ring |6(] and magnetic molecular clusters in 
longitudinal and transverse ac magnetic fields [p.9| . 

We consider the cooperative effects of two monochro- 
matic fields modulating the bias energy he and the TME 
hA. The dissipative TSS is described by the spin-boson 
Hamiltonian p|,^,pT[ (we set H = fc B = 1) 



H(t) = -i[ e (t) az + A(t)a x ] 
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where Hb = J2i[Pi/^ rn i + rn i UJ i x i I "A represents the bath 
of bosons, and the collective variable X = J^\ CiXi de- 
scribes the bath polarization. All effects of the boson 
bath on the TSS are captured by the spectral density 

2 

J(lu) = ? J^. — i — S(u> — ojj). A dipole interaction contri- 
butes to the bias energy, whereas a quadrupole modula- 
tion results in an exponential modification of the TME, 

A(i) = A exp[A A sin(fi A t)] , e(t) = e 1 sin(fi e i) . (2) 

In the absence of driving, A& — 0, e 1 = 0, the equilib- 
rium state reached at long times has the localized states 
\R) and \L) occupied with equal probability. If either the 
bias or the coupling energy is modulated, the left-right 
symmetry is dynamically broken. However, on average 
over a period, one finds asymptotically again equal occu- 
pation of both states. When both parameters are modu- 
lated with commensurable frequencies (n, m integer) 



a 



mfl and fL 
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the Hamiltonian p) is periodic (we assume T = 2n/fl to 
be the smallest period in common). We find that equal 
occupation on average is reached at long times when n 
or m is even. However, when both n and m are odd, 
the right-left symmetry is broken even on average. The 
symmetry breaking is maximal when n — m = 1. The 
phenomenon that the occupation of one of the two states 
is preferred on average at long times is schematically 
sketched and explained in Fig. [y. In the following, we 
give a quantitative study of this effect. 

Suppose that the TSS has been prepared at time zero 
in the state \R) (<j z = +1) with the bath in thermal 
equilibrium. The dynamical quantity of interest is then 
the population difference (cr z ) t = P{i) = PR,{t) — Pl(^) 
for this factorizing initial condition. Since the bath is 
harmonic, it can be traced out exactly. Then, P(t) is 
expressed as a double path integral over the forward and 
backward spin paths <j(t) and <j'{t) which are piecewise 
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t = 3774 




FIG. 1. The bistable potential (top) and the related TSS 
(bottom) are symbolically sketched at different times for the 
case Q = Q A = fi e = 2-k/T. At time t = 0, the TSS is 
symmetric. At time t = T/4, the right state is lower and the 
TME has the maximum value. At time t = ST/4, the left 
state is lower, but the TME has the minimum value. Thus, 
relaxation towards the right state is preferred. 

constant with values ±1. The effects of the environment 
are in an influence functional which introduces nonlocal 
correlations between different path segments g. It is 
convenient to switch to the path combinations 77 (r) = 
[ct(t) + o-'(r)]/2 and £(r) = [a(r) - ct'(t)]/2. Then, P(t) 
is expressed in terms of the double path sum 
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where A is the path weight in the absence of the bath 
coupling, and the influence function is 
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FIG. 2. Plots of the population difference P(t) for a — 1/2 
as calculated from the expression (0). The driving-induced 
symmetry breaking phenomenon is clearly visible. The pa- 
rameters are T = 0.05A , uj c — 50A , e 1 = 5A , and 
Q t = S1 A = 5A . The upper curve describes the case 
Aa — 0.5, and the lower curve the case Aa = —0.5. The 
horizontal lines follow from Eq. ([13) with (FLjEJ) . The middle 
curve corresponds to a time-independent TME (Aa — 0). In 
this case, the right/left symmetry is asymptotically restored. 

For the special case a — 1/2 of the Ohmic coupling, 
however, the dynamics can be found in analytic form in 
the scaling limit (pj) . Using the concept of collapsed blips 
@, we find 



x [Q'(h - tjfa) + iQ"(h - hWh)] . (4) P s (t) = e~ G ^ , G(t 2:tl ) 
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dsA 2 (s) , 



For the Ohmic spectral density J (to) — 2itauj e~ u / Uc , 
where a is the dimensionless Ohmic coupling, the kernel 
takes for times lu c t 3> 1 the so-called scaling form 

Q(t) = 2a\n[((3u! c /n) sinh (7rr//?)] + i7rasgnT . (5) 

It is convenient to write P(t) = P s (t) + P a (t), where 
P s (t) and P a (t) are the symmetric and antisymmet- 
ric parts under inversion of the bias e —> —e. Since 
limj^^ P s (t) ~ 0, the asymptotic value is determined 
by the antisymmetric part P a (t). In the absence of the 
biasing force, t x = 0, we have P a {t) = for all t. For any 
initial preparation, we then have \ira t ^ oa P(t) = 0, i.e., 
the TSS is completely delocalized at long times. Simi- 
larly, if we choose Aa = and e x ^ 0, the time average 
(sin('(i 2 ,i 1 )) r over a period 2n/Q, e of the external field 
vanishes. Thus, on average, the TSS is again delocalized, 
i.e., P w = lim t _ +00 (P(i)) r = 0. However, the construc- 
tive interference between diagonal and off-diagonal driv- 
ing, i.e. A/a ^ and e\ ^ 0, can lead to driving-induced 
symmetry breaking, as qualitatively illustrated in Fig. \U. 
We find P m ^ for f2 e = fi A , and the tight-binding site 
preferred is determined by the relative sign of A/\ and t x . 

The evolution of P(t) can not be given in analytic form 
apart from special limits. Numerical methods cover short 
to intermediate times, but they become costly or even in- 
adequate if one is interested in the asymptotic behavior. 



P a (t)= / dt 2 / dt x e 
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x e- Q '^-^ sin [dt^t,)] e - G (*2>*i)/2 



with C,(t 2 , £1) = e-L f. 2 ds sm(tt e s) . The relaxation of P(t) 
is shown for different A A in Fig. 2. The superimposed os- 
cillation is due to the oscillatory bias. The middle curve is 
for A A = 0. In this case, the TSS is asymptotically delo- 
calized on average. Symmetry breaking only occurs when 
both the bias and the intersite coupling are monochro- 
matically modulated (upper and lower curve). The re- 
sults depicted in Fig. 2 correspond to the case fi A = fl e . 
Evidently, the sign of P a (t) depends on the sign of A/\. 

In the general case, the dynamics of the driven system 
can be described by the exact master equation |20) 

P(t)= f dt'[K a {t,t')~K s (t,t')P(t')}. (7) 

Jo 

The familiar noninteracting-blip approximation (NIB A) 
describes the dynamics correctly for high temperature 
and/or strong damping J3J3]. The NIBA kernels read 



K s/a (t,t 



'\ - P -Q'(t-t': 



B s/a [Q"(t-t')]L s/a (t,t'), (8) 



where B s (z) = cos(z) and B a (z) = sin(z). The effects of 
the driving forces are captured by 



L s/a (t,t<) = A(t)A(t')B s/a [at,t')} 



(9) 



The equilibrium state for a TSS with static bias e is 



[k + (e) - fc_(e)]/[fc + (e) + fc -(e)L where fc ±( e ) 



IS 



the forward/backward rate for the static case. One has 
k±(e) = (Ag/4) j^ dt exp{±iet—Q(t)], obeying detailed 
balance, k_(e) = cxp(— j3e)k Jt 
Ohmic forward rate reads M 



(e). In the scaling limit, the 



k+{e) -4u Jc \2n) r(2a) ' ( j 

For the driven TSS, Eqs. (|l|) - (||), the average rate is 



k± = 



2^ 



1 



2-rr/n 

dtl dT^[K t ll.l-,) = K„(l.l-,) 



Next, we expand the bias phase terms e T*(«i/^«) cos ( fi «*) 
and e ± J ( e i/ n O cos [ n =( f — n] m ser i es i n J-Bessel functions, 
and the oscillatory factors of the tunneling coupling 

e A A sin(fi A t) and e A A sin[n A (t-r)] mto g^g m /_B eS Sel 

functions. We then find the average rate k± in the form 




FIG. 3. The mean asymptotic value Poo is shown as 
a function of the bias amplitude e x for Ohmic couplings 
a = 0.1 (full curve), a — 1/4 (long-dashed curve), and 
a = 1/2 (short-dashed curve). The parameters are Aa = 0.5, 
il c = fi A = A /2, T = A , and oj c = 50A . The inset 
shows Pqo versus a for temperatures T = 0.1 A (short-dashed 
curve), T = A (dashed curve), and T = 10 A (full curve). 
The bias amplitude is e 1 — 10 A . 



g e fcr(n/m-l)(r+.)/2 fc ± [( m ^_ ns )ft] 



X / r (AA)/ s (A A )J^ — J Jt-( r +s)n/m 
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(11) 



The individual channels in (|llj) describe tunneling under 
emission and absorption of a fixed number of quanta with 
frequencies which are multiples of the frequency Q. The 
channel weights are given in terms of / and J Bessel 
functions. The exclusive channel rates are expressed in 
terms of the rates fe ± for a static bias e = (ml — ns)fl. 
The mean asymptotic position is then found as 



P„ 



(k + -k_)/(k + + k_). 



(12) 



Substituting Eq. ( [JT|) and the symmetry relation fc_(e) = 
fc+(— e), and using the properties J_„(z) = (— l)™J n (z) 
and I- n (z) = I n (z), we find P^ = when n or m is 
even, whereas Poo =/= when both n and m are odd. 

The predictions following from ( |l2| ) are shown in Fig. g 
(horizontal lines) in comparison with the exact P(t) (os- 
cillating curves). We have chosen a — 1/2 to allow for a 
direct comparison with the predictions of the exact solu- 
tion (ra). Evidently, the asymptotic position value is very 
accurately reproduced by the NIBA. 

In Fig. pi we show P m as a function of the bias am- 
plitude e x for different coupling strengths. The full curve 
(a = 0.1) exhibits oscillations. The oscillations are al- 
ready damped out in the curves for the higher damping 
values displayed. Interestingly, the symmetry breaking 
effect grows with increasing a and increasing e x . For 
large bias amplitude, in the investigated bias regime, P^ 
saturates. The inset depicts P M versus a for e 1 = 10A 
and different temperatures. For weak coupling a C 1, 
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FIG. 4. Plots of Poo versus bias e 1 for fi A = Q t — A /2 
(full curve), fi A = 3f2 e = A /2 (long-dashed curve), and 
fi,, = 3fi A = A /2 (short-dashed curve). The other parame- 
ters are a = 1/4, Aa = 0.5, T = A and uj c — 50A . 

the symmetry breaking becomes more effective with in- 
creasing temperature. For stronger coupling, on the con- 
trary, it is more pronounced the lower the temperature. 

In Fig. EL we show P M for different ratios fl € /H,A- The 
full and long-dashed curves represent the cases fl e = fl A 
and £1a = 30 e - In both cases, the symmetry breaking is 
clearly evident. When the bias frequency is three times 
the barrier modulation frequency, the symmetry break- 
ing effect is strongly reduced (short-dashed curve). The 
different behavior is due to the oscillatory and smooth 
behavior of the J n (z) and I n (z) functions, respectively. 

We believe that the phenomenon of dynamical symme- 
try breaking can be seen, e.g, in Josephson systems. A 
promising candidate is a superconducting quantum inter- 
ference device in which the tunneling coupling and bias 
energy can be varied independently H. In the variable 



barrier rf SQUID, the single Josephson junction in a stan- 
dard rf SQUID ring with inductance L is replaced by a 
superconducting loop with inductance Ld c interrupted by 
two junctions. In the limit Ld c /L — > 0, the 2D Joseph- 
son potential V($,&dc) goes over into the ID potential 
of a standard rf SQUID with adjustable critical current 
Ic($xdc)- Within the RSJ model, the deterministic equa- 
tion of motion of the rf SQUID for the flux variable 
<t> = 2tt$/$ = 2e$ is Cij> + j>/R + 4e 2 (d / 'd(f>)U (</>) = 0. 
The rf SQUID potential is given by {U = l/4e 2 L) 



U(4>) = U [ 



0x) 2 -/3L(^dc)cos(0)], (13) 



where ^ti^xdc) — ^ e LI c cos(^c/) x dc)- For an applied flux 
4> x = 7r and /3l > 1, the potential is symmetric about 
<j> = 7r with two degenerate minima at <j> = 7r ± m , where 
m is the smallest positive root of <f> m — /3_Lsin0 m = 0. 
The barrier height is Ub = t/o [/3l — ^m ~ Pl cos </> m ] , and 
the squared well frequency is u)q = (1 — /?l cos <f> m ) / LC . 

With the modulation 4> x = it + 1 sin(f7 e i), where 
<t>i *C 0m, the system is effectively described at low T by 
Eq. (|lj) with the bias variation (|2|), and amplitude e 1 = 
(j) m (j) 1 /2e 2 L. The parameter w = Ub/uio can be varied m- 
dependently of the bias by varying the flux <p x d c applied 
to the small loop. For (3l ^> 1, the tunneling length 20 m , 
and hence the damping parameter a = (<fi m / ir) 2 Rq / R 
(Rq = 7r/2e 2 = 6.5 kCl), is approximately constant un- 
der variation of <p x dc- Harmonic modulation of the WKB 
exponent v of the TME A oc exp(— cv), where c depends 
on the barrier shape, can now be achieved by driving 
(frxdc- Thus one can experimentally reach the situation 
described by Eq. (ph. Typical realizable junction and 
loop parameters are L w 200 pH, C rs 100 fF, /3l ~ 5, 
and R/Rq ranging from 0.5 to 10. This results in values 
for a ranging from 0.1 to 2, barrier height Ub ~ 1 ... 10 K, 
well frequency cu Q w 10 10 . . . 10 11 s" 1 , and the tunneling 
matrix element in the MHz regime. We expect that the 
symmetry breaking effect can be observed at tempera- 
tures about 10 mK and driving frequencies ^a and Q e 
in the 1 ... 10 MHz regime. 

A giant spin or nanomagnet is another possible candi- 
date for the observation of the above symmetry breaking 
phenomenon. In these systems, the bias can be tuned by 
a magnetic field along the easy axis. It has been demon- 
strated that the TME of the Fe 8 molecular system can be 
tuned by a transverse magnetic field |H| , as theoretically 
predicted p2| . Therefore, it seems feasible to operate 
nanomagnets in the appropriate parameter regime. 

In summary, we have found that dynamical symmetry 
breaking can be induced by a concerted modulation of the 
bias and tunneling energy of a symmetric TSS. Because of 
constructive interference, dwelling in one of the localized 
states is preferred on average. This effect is robust since 
it is largely independent of the dissipative mechanism. 
Also, the dependence on coupling strength and temper- 
ature is interesting. For weak coupling, the asymptotic 



mean position increases with temperature, whereas for 
strong coupling it shows opposite behavior. 

Readily, the results can be extended to quantum trans- 
port in a tight-binding (TB) lattice. In the incoherent 
tunneling regime, the asymptotic current is proportional 
to the asymptotic mean position of the TB system [pTl . 
Upon imposing concerted driving of bias and TME, it is 
possible to extract a finite current out of an (on average) 
undirectional force. Thus, the system constitutes a novel 
realization of a discrete quantum Brownian rectifier pl| . 
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